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Abstract. A new class of simple and exact solutions of relativistic hydrodynamics 
is presented, and the consequences are explored in data analysis. The effects of 
longitudinal work and acceleration are taken into account in an advanced estimate 
of the initial energy density and the life-time of the reaction. 



1. Introduction: High temperature superfluidity 

Fluid dynamics is based only on local conservation laws and thermodynamics, so it can 
be successfully applied to a vast range of physical phenomena. Recently, a new state of 
matter has been created in Au-|-Au collisions at RHIC, and, surprisingly, it was found 
to flow as a perfect fluid [1]. The kinematic shear viscosity of this nearly perfect fluid 
has been determined and found to be at least a factor of 4 smaller, than that of the 
superfluid ^He [2]. Let us refer to this property as high temperature superfluidity, to 
highlight that the fluid of quarks in the super-high, T ~ 2 terakelvin temperature range 
flows better than the superfluid ^He in the extremely low, 1-4 K temperature region. 

Exact solutions of perfect fluid hydrodynamics have powerful implications in high 
energy particle and nuclear physics. Although the renowned Landau-Khalatnikov 
solution [3, 4, 5] is only implicitly given, it describes well the energy dependent increase of 
the width of the rapidity distribution. Although the renowned Hwa-Bjorken solution [6] 
lacks acceleration and yields a too perfectly flat rapidity distribution, it yields a key 
estimate of the initial energy density in high energy reactions. Other families of exact 
solutions were born from the desire of understanding the dynamics of high-energy heavy 
reactions [7, 8, 9, 10, 11, 12, 13, 14]. Here we present a new, accelerating, analytic, exact 
and explicit solution of relativistic hydrodynamics, which flts dn/dy data at RHIC and 
allows for an advanced estimate of initial energy density and life-time of the reaction. 



2. Simple solutions of relativistic hydrodynamics 

We discuss spherical (and one-dimensional) solutions, where r denotes the radial (or 
the single spatial) coordinate, and d is the number of dimensions. The pressure is 
denoted by p , the energy density by e , and the temperature by T, the four-velocity 
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Table 1. New exact hydrodynamical solutions arc given by lines (a, b, c) and (e), 
while case (d) is the already known Hwa-Bjorken (and Hubble) solution. 



by u*^ — 7(1, v), and v = vn is the three- velocity. For the equation of state (EoS), 
we assume that e — B — k{p + B), with k — l/c^, where Cg is the speed of sound and 
B stands for a bag constant (that may have a vanishing value too). In high energy 
collisions, the entropy density a is large, but the net charge density is small. Thus we 
assume that all the conserved charges have zero chemical potential. In perfect fluids, a 
and four-momentum tensor T^^ = {e + p)u'^Ui, —p5^ are locally conserved: du{au^) = 0, 
dyT^^ = 0. With projections, we obtain the Euler and the energy equations: 

{s + pXd.u'' = {S^ - u'^Up) d^p, (1) 

{e + p)d,u'' + u''d,e^O. (2) 

Our solutions can be written down in the so-called Rindler-coordinates r and 77: 
t — T cosh 77, r = TsinhTy. We found the following expression for p and v: 

V = tanh Xr] , p + B = po ^ — j " ^cosh - j . (3) 

We obtain solutions in some special cases of the constants listed in Table 1, and for any 
B. In what follows, we set -B = to simplify the description. While case (d) is the 
Hwa-Bjorken solution, case (a) is a new solution valid in arbitrary dimensions. It has 
uniformly accelerating trajectories. Case (b) and (c) were found jointly by T. S. Biro [15] 
(for d — 1) and by us (for any d). In these solutions p is finite in 77. Case (e) has a 
special EoS and is valid only in 1-1-1 dimensions, but A can be any real number. This is 
why we shall use it in the following applications. 

3. Applications to high-energy reactions 

For the description of high-energy reactions, we calculated the rapidity distribution, 
^ from the solution denoted as case (e) in Table I. Assuming sudden freeze-out at 
temperature T{Tf,r] — 0) — Tf on a, hypersurface pseudo-orthogonal to u^{x), with 
a saddle-point integration in 77 and the transverse mass rriT (which becomes exact if 
m/Tf » 1, where m is the particle mass), we got the following formula [20]: 

cosh-t- m exp 1-^ fcosh" m - ll 1 , (4) 
dy dy y=o \aJ ^ \ Tf I \aJ Jj' 
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Figure 1. Left panel: Result of the fit of eq. (4) to BRAHMS ^ data of Ref. [16]. 
Right panel: Flow trajectories and possible freeze-out hypersurfaces of the A ~ 1.18 
solution, which fits to BRAHMS ^ data. (The physical solution is in the future light 
cone, but for aesthetic reasons its extension to the past light cone is also shown). 



with a = ^rj-- The parameter Ay"^ = ^jT^jXl^+iia characterizes this distribution: 
it has a minimum at ?/ = 0, if Ay^ < 0, it is flat if Ay^ = 0, (this is the case 
when A = 1), otherwise it is nearly Gaussian. We extracted the A parameter for 
collisions at ^/sNN = 200 GeV by fitting eq. (4) to g data measured by the BRAHMS 
collaboration [16]. Fig. 1 shows the fit, and the corresponding flow profile. We found 
A = 1.18 ± 0.01, which indicates the presence of acceleration, as shown in Fig 1 by the 
curvature of the fluid lines. This influences both the estimation of initial energy density 
(because of the faster expansion and the work done by the matter), and the estimation 
of the life-time of the reaction. This way we found a new estimate of the initial energy 
density Ec that generalizes Bjorken's estimation, ebj to non-flat rapidity distributions: 

^ = (2A-l)(2)"' where ^b, = . (5) 

EBj \TqJ {R^7l)Tody y=0 

The life-time estimation from the longitudinal length of homogeneity, -Riong of Sinyukov 
and Makhlin, tbj was also based on the Hwa-Bjorken solution [18]. Our new estimation, 
Tc, which takes acceleration effects into account, is 

-Rlong = a/— T = Atbj = A./^i?long. (6) 

y nit X V ^/ 

The eqs. (5) and (6) show that both the initial energy density and the life-time of the 
reaction is under-estimated by formulas based on the accelerationless of Hwa-Bjorken 
solution. For more physical EoS than that super-hard k = 1 corresponding to the 
A = 1.18 solution, we conjectured new, advanced estimations of the initial energy density 
Er2 and the life-time of the reaction as 

'-S 

i^i=(2A-l)(lt) m , (7) 

=[X+(X-l)(l-cl)]TBr (8) 

For their detailed explanation, see Refs. [19, 20]. Using the Bjorken estimate of Ebj =5 
GeV/fm^ as given in Ref. [17], and tj/tq = 8 ± 2 fm/c, our analytic formula (5) yields 
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Figure 2. Dependence of Ec^/eBj on t//to (left panel) and Tc2/TBj on A (right panel) 
for various values. Solid lines for — 1 indicate the exact result of cqs. (5) and (6), 
the others show the conjectures of eq. (7). 



Ec = (2.0 ± 0.1)eBj = 10.0 ± 0.5 GeV/fm^ for A = 1.18 ± 0.01. For the life-time, eq. (6) 
implies a 18±1% increase. Fig. 2 shows the speed of sound dependence of Sci/sBj as 
a function of tj/tq, and that of r^^jTsj as a function of A. For a reahstic equation of 
state, where = 0.1 in ^Jsmi = 200 GeV Au+Au reactions [23], acceleration increases 
the initial energy density estimate to 15 GeV/fm^, and the life-time estimate by 36 %. 
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